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ABSTRACT 

The comparison between the new Spitzer data at 24 /im and the previous ISOCAM data at 15 /jm is 
a key tool to understand galaxy properties and evolution in the infrared and to interpret the observed 
number counts, since the combination of Spitzer with the Infrared Space Observatory cosmological 
surveys provides for the first time the direct view of the universe in the infrared up to z>2. We 
present the prediction in the Spitzer 24 /im. band of a phenomenological model for galaxy evolution 
derived from the 15 /im data. Without any "a posteriori" update, the model predictions seem to agree 
well with the recently published 24 /im extragalactic source counts, suggesting that the peak in the 
24 /im counts is dominated by starburst galaxies like those detected by ISOCAM at 15 /im but at 
higher redshifts (1 1 z t 2 instead of 0.5 < z < 1.5). 

Subject headings: galaxies: evolution - galaxies: starburst - cosmology: observations - infrared: 
galaxies. 



1. INTRODUCTION 

Cosmological constraints on the evolution of galaxies 
have been recently investigated by studying the statisti- 
cal properties of large samples. In particular, the com- 
bined analysis of extragalactic source counts and redshift 
distributions is generally used to calibrate the theoretical 
models as a function of time. The mid-infrared (MIR) 
and far-infrared (FIR) regions of the electromagnetic 
spectrum probe the population of actively star-forming 
galaxies obscured by dust. Extragalactic source counts 
from different surveys over a wide flux range obtained 
with the ISOCAM instrument on board of the Infrared 
Space Observatory (ISO) indicate that these sources have 
evolved rapidly, significantly faster than deduced from 
optical surveys (Elbaz et al. 1999; Gruppioni et al. 2002; 
Metcalfe et al. 2003; Rodighiero et al. 2004). These re- 
sults are supported by the detection of a substantial cos- 
mic infrared background (CIRB; Hauser & Dwek 2001), 
which is interpreted as the integrated emission from dust 
present in galaxies. The contribution of resolved ISO- 
CAM sources accounts for -60-70% of the CIRB at MIR 
frequencies, the bulk of this background originating in 
discrete sources at z<1.2 (Franceschini et al. 2001; El- 
baz et al. 2002). 

The Spitzer Space Telescope is now providing new in- 
sights into the IR population contributing to the CIRB, 
in particular with the Multiband Imaging Photometer 
24 /im band, which is starting to detect a population of 
galaxies that may be IR-luminous galaxies at z — 1.5 — 3 
(i.e. high-redshift analogs of the faint 15 /im galaxies 
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detected by ISOCAM; see Chary et al. 2004). Prelimi- 
nary results from the 24 /im extragalactic source counts 
(Marleau et al. 2004; Papovich et al. 2004), confirm- 
ing the existence of the rapidly evolving dust-obscured 
population discovered bylSOCAM, raise the question on 
how to compare them with the previous ISOCAM counts 
at 15 /tm. Both 24 and 15 /im bands are extremely sen- 
sitive to the presence of broad emission features at 6.2, 
7.7, 8.6, 11.3 and 12.7 /im in the spectra of galaxies, 
probably from polycyclic aromatic hydrocarbons (PAHs; 
Puget & Leger 1989). Since these features dominate 
the photometric output at some redshifts, the ratio be- 
tween the Spitzer 24 /im and the ISOCAM 15 /im fluxes 
(<S , 24/j.m/ 'Sisum) is strongly dependent on z. In figure Q 
this ratio versus z is shown for the populations contribut- 
ing to the MIR source counts (see next section for a dis- 
cussion): starburst galaxies (M82 spectral energy dis- 
tribution [SED]), normal galaxies (M51), type 1 active 
galactic nuclei (AGNs; SED from Elvis et al. 1994), and 
type 2 AGNs (Circinus). As clearly shown in the plot, the 
15 /im band is optimized for detecting 0.5 <z< 1.5 galax- 
ies, while the 24 /im band is favored for the detection 
of galaxies at z> 1.5. Since the comparison between the 
source counts in the two bands is a powerful tool for un- 
derstanding the evolutionary properties of the different 
galaxy populations contributing to the counts at differ- 
ent redshifts and flux levels, it is worthwhile performing 
a careful comparison. 

In a recent paper on extragalactic source counts from 
the First Look Survey (FLS), Marleau et al. (2004) try 
to compare the 24 /im source counts with the previous 
15 jum counts from different ISOCAM surveys, trans- 
formed to 24 /im. However, the reported transformation 
appears to be overly simplistic, since the ISOCAM counts 
(plotted in figure 4 of Marleau et al. (2004)) have been 
converted to the Spitzer 24 /im band by considering a 
single value for the SW/xm/ 'S^m flux ratio for all flux 
densities. The value used by the authors (Ss^m/ Sis^m 
= 1.2) is a median value derived from typical luminous 
infrared galaxy/ultraluminous infrared galaxy SEDs at 
relatively high z (Chary & Elbaz 2001), which is ap- 
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propriate only for sources making up the peak of the 
counts. In fact, as shown in figure ^ if we consider the 
"starburst" template (M82), we obtain a local value for 
the flux ratio of ~2.5, while we reach a value of ~1.3 
only at z ~ 1. Therefore, only the contribution to the 
counts made by galaxies with z around 1 could be trans- 
formed to 24 /im using a flux ratio similar to that consid- 
ered by Marleau et al. (2004). In particular, the bright 
part of the European Large-Area ISO Survey (ELAIS) 
counts (Gruppioni et al. 2002), linking the IRAS counts 
to the deep ISOCAM counts, is dominated by nearby 
non-evolving normal galaxies (with ratios of ~1.7) and 
by starburst galaxies and type 2 AGNs (with ratios of 
~2.5 and ~2.3, respectively), as shown by La Franca et 
al. (2004). Therefore, the use of a single ratio value of 
1.2 produces a misleading result [especially at 5 , 24 Aim > 1 
mjy: i.e. 15 /im counts shifted downwards by a factor 
of (1.7/1.2) 15 - (2.5/1. 2) 15 = 1.7 - 3.0], suggesting an 
apparent inconsistency between the bright part of the 
source counts in the two bands. 

Since at the moment there are no large areas with avail- 
able data at both wavelengths, we can make use of a 
model fitting the observed 15 /jm source counts (Pozzi 
et al. 2004, I. Matute et al. 2005, in preparation) to 
transform the counts from one frequency to the other, 
thus allowing a direct comparison. 

In this Letter we discuss the more realistic way to 
transform the model predictions and the observed data 
from the LW3 band of ISOCAM to the 24 /im band of 
Spitzer, in order to compare the properties of the 24 
/im sources with those of the ISOCAM 15 /im ones. The 
Letter is structured as follows: in sectional we describe 
the evolutionary model fitting the 15 /im observables; in 
section l3~Tl we show the model predictions at 24 /im; in 
section 13.21 we discuss a method to transform the ob- 
served data points from 15 to 24 /im; in section 0] we 
present our conclusions. 

Throughout this Letter we will assume H$ = 75 km 
s" 1 Mpc"\ n m = 0.3 and fl A = 0.7. 

2. THE MODEL 

The model is based on the first direct determination 
of the 15 /zm luminosity function (LF) of galaxies and 
AGNs, based on data from the ELAIS southern fields 
survey (Lari et al. 2001; La Franca et al. 2004, Rowan- 
Robinson et al. 2004). We assume that four main pop- 
ulations, evolving independently, contribute to the ob- 
served source counts: starburst and normal galaxies and 
type 1 and 2 AGN. A maximum likelihood analysis (Mar- 
shall et al. 1983) has been used to simultaneously fit 
both evolution rates and shape parameters of the differ- 
ent local LFs. Although AGN make up a non negligible 
fraction of the extragalactic source counts at 15 /im (es- 
pecially at high flux densities), galaxies are the dominant 
class in the MIR. 

The LF determination for galaxies, described exten- 
sively in Pozzi et al. (2004), is based on a sample of 
about 150 galaxies in the redshift interval 0.0 < z < 0.4, 
covering a large flux density range between I/MS and the 
deep ISOCAM surveys (0.5 < S 15tim < 50 mjy). The 
normal, non evolving, and the starburst, evolving, pop- 
ulations are separated using the new criterion based on 
the MIR to optical luminosity ratio (Iis^m/Ifi)- We use 
the basic Silva et al. (1998) models to reproduce the SED 



of our prototypical galaxies, assumed to be M82 for the 
starburst population and M51 for the normal one. The 
MIR region (between 3 and 18 /im) of the modeled spec- 
tra were replaced with ISOCAM circular variable filter 
observations (M82: Forster-Schreiber et al. 2003; M51: 
Roussel et al. 2001). Note that, for simplicity, we have 
used a single template SED for each population, instead 
of different SEDs for different infrared luminosity inter- 
vals. While the normal population is consistent with no 
evolution, for the starburst population a strong evolution 
is found both in luminosity [L(z) oc (1 + z) 3 ' 5 up to z ~ 1] 
and in density [p{z) oc (1 + z) 3 8 up to z ~ 1], The evo- 
lutionary parameters of our model have been tested by 
comparing the model predictions with all the available 
observables, like source counts at all flux density levels 
(from 0.1 to 300 mjy) and redshift distributions and LF 
at high z. The agreement between the model predictions 
and the observed data is remarkably good (see figure [5] 
for an example of how the model fits the observed num- 
ber counts at 15 /im ). 

The LF determination for AGNs (both type 1 and 2), 
described in I. Matute et al. (2005, in preparation), is 
based on ELAIS data (27 type 1 AGNs and 25 type 2 
AGNs) combined with the local IRAS sample at 12 /im 
of Rush, Malkan & Spinoglio (1993), converted to 15 /im 
using appropriate SEDs (41 type 1 AGNs and 50 type 
2 AGNs). The typical SED assumed for type 1 AGNs 
is that compiled by Elvis et al. (1994), while for type 
2 AGN two extreme cases of the obscured AGN SED in 
the MIR were assumed: a starburst-likc SED (Circinus; 
Sturm et al. 2000) and an AGN-like SED (NGC 1068; 
Sturm et al. 2000). Type 1 AGNs are found to evolve 
with a luminosity evolution [L(z) oc (1 + z) L )], with 
an evolution rate kh equal to 2.6 up to z ~ 2 and con- 
stant thereon. A similar evolutionary scenario is found 
for type 2 AGNs, with Ul ranging from 2.0 to 2.6 depend- 
ing on the adopted SED (NGC 1068 or Circinus, respec- 
tively) . The best-fitting model is found to reproduce well 
both observed source counts and redshift distributions, 
as shown by I. Matute et al. (2005, in preparation). 

3. COUNTS AT 24 MICRON 

3.1. Model Predictions at 2^ fim 

By using the appropriate SED for each population and 
convolving the SED with the appropriate filter transmis- 
sion, the galaxy and AGN local LFs have been trans- 
formed from 15 to 24 /im. Then the predicted 24 /im 
source counts have been computed for all the contribut- 
ing populations. In figure |21 the source counts pre- 
dicted by our model are compared to the recently pub- 
lished Spitzer 24 /im data from the FLS (Marleau et al. 
2004), from the deep surveys (Papovich et al. 2004) and 
from the Spitzer Wide- Area Infra-red Extragalactic Sur- 
vey (SWIRE; D.L. Shupe et al. 2005, in preparation). 
We can notice a consistency between data and model, 
with no need, at least in first approximation, for the use 
of more extreme SEDs for starburst galaxies (i.e., Arp 
220). AGNs (either type 1 or 2) do not dominate the 
observed source counts at any flux level, although type 2 
make about 25% of the counts at SW/xm > 10 mjy. The 
counts arc dominated by non evolving normal galaxies at 
> 8 mjy and by evolving starburst galaxies at lower flux 
densities. The assumption made in our model of no evo- 
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lution for galaxies at z > 1 , not very well constrained by 
ISOCAM data, is "a posteriori" consistent with Spitzer 
data. 

It is interesting to see how the ratio between the Spitzer 
24 and the ISOCAM 15 fim flux for all the popula- 
tions changes as a function of z (figure ^1 and of the 
24 fim flux (as derived by our model; figure HJ . The 
comparison between figs ^ and 0] clearly shows that the 
higher flux densities (S^^m > 2 — 3 mJy) are dom- 
inated by nearby objects with moderately high values 
of the 5 , 24 A im/5'i5/im ratio (~ 2.2: starburst and type 2 
AGNs; ~ 1.7: normal galaxies), while the bump of the 
24 fim counts (at fluxes 0.1-1 mJy) is dominated by ob- 
jects with 5 , 24 Aim /S'i5 Atm « 1.4 (mainly starburst galaxies 
at 0.7 < z < 1.5). These are the same populations found to 
contribute to the ISOCAM 15 /urn source counts. How- 
ever, since we have approximately 
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and the ratio value at flux densities around the peak of 
the differential source counts is ~ 1.4, the evolutionary 
excess is more pronounced at 15 /im than at 24 /im. At 
the lower flux densities (SW^n ~ 0.1 mJy), sources with 
high flux ratio values (5 , 24 AI m/S'i5 Aim > 2.0 — 2.5, cor- 
responding to z > 1.5) start to dominate the counts. 
These high-z starburst galaxies, mainly at 1.5 < z < 2.5, 
are not visible in the ISOCAM surveys, because of the 
15 fim fc-correction, but contribute to the fainter end of 
the 24 fim source counts and to the 24 /im cosmic back- 
ground. Since the existence of this starburst population 
at high redshift is predicted just by extrapolating our 15 
/im model to higher z, it is likely that faint 24 /jm sources 
are the high- redshift likes of the fainter 15 fim sources de- 
tected by ISOCAM (see also Chary et al. 2004). 
The same kind of considerations are evident from figure 
where the contributions to the 24 /im number counts 
from different redshift intervals are shown. While galax- 
ies with 0.0 ^z;; 0.5 dominate the high fluxes (S24 Mm > 2 
mJy), the peak of the differential counts is made manly 
by 0.5^2:^1.5 sources, although at 5 , 24 /jm J: 0.16 mJy 
(corresponding roughly to the limit of the deepest, non- 
lensed, ISOCAM survey: Si 5fJim = 0.12 mJy; see figure 
|SJ) the higher- z population contributes for a conspicuous 
fraction (~46%). By extrapolating the source counts 
model down to the faintest fluxes, we obtain an esti- 
mate of the total 24 /im background: v I „(2A fim) ^2 A 
nWm _2 sr _1 . According to our model, the amount pro- 
duced by sources at z< 1.5 is ~1.6 nWm~ 2 sr _1 , therefore 
^66% of the 24 /im background originates at relatively 
low redshift. From the comparison between the total 
background predicted by our model and the value de- 
rived from the observed 24 /im source counts (1.9 ± 0.6 
nW m~ 2 sr^ 1 ; Papovich et al. 2004), we find that the 
Spitzer deep surveys already resolve ~78% of the total 
24 /im background. 

3.2. Observed ISOCAM Counts Transformed to 24 fim 

To compare observed data counts at 15 fim with 
those at 24 fim, we have converted the 15 /im source 
counts obtained from different ISOCAM surveys (from 
the ultra-deep lensed of Metcalfe et al. 2003, to the shal- 
lower ELAIS Survey of Gruppioni et al. 2002) to 24 



/im, as described in the following text. We have con- 
volved the observed 15 /im differential counts (- 



with the distribution of the ratios S^^m/Sis^m ob- 
tained from our model, given a 15 /im source selection 

(/ [(£24^771/ Sl5/JTTi), Sis^m]): 
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Data at 15 /im from different samples have been com- 
bined by weighting each point by its formal error (inverse 
of the squared error). In figure|21the counts derived from 
the 15 /im observed data (shaded area) are overplotted 
to the 24 /im data and model. The two different source 
counts are consistent and, in first approximation, seem 
to agree well. In particular, at high flux densities we do 
not observe the discrepancy as reported by Marleau et 
al. (2004), thanks to a more accurate flux density ratio 
applied. On the other hand, the discrepancy observed 
at 24 fim fluxes lower than ~0.05 mJy is only apparent, 
since fluxes fainter than this are not sampled by the 15 
/im data. Some level of inconsistency between the two 
source counts are visible around ^1 mJy, where the 15 
fim counts (and the model) are higher than the observed 
24 fim source counts. A decrease of the S^/jm/SiS/im ra ~ 
tio around 1 mJy could be obtained by slightly modifying 
the starburst template in the MIR range (i.e. increasing 
the PAH features with respect to the continuum, since 
at the typical redshift of sources with these flux densi- 
ties, z ~ 1, the PAH features enter the 15 fim band). 
A similar change was made by Lagache et al. (2004), 
who show how a minor change in the starburst template 
spectra between 12 and 30 fim (together with a slight 
modification of the luminosity density) was a sufficient a 
posteriori modification to enable a model not fitting the 
observed 24 fim source counts (Lagache et al. 2003) to 
reproduce the observations. 

We are actually working at improving the model-data 
agreement by including also the recently published 24 fim 
source counts (and all the MIR-FIR counts and redshift 
distributions available in literature) as an observational 
constraint to the maximum likelihood analysis of our 15 
/im-based model. Moreover, we are considering the use of 
an SED library, with different SEDs associated not only 
to different populations, but also to different luminosity 
classes (i.e. Chary & Elbaz 2001). The model improve- 
ments will be described in a future paper (F. Pozzi et al. 
2005, in preparation), since they are beyond the scope of 
the present Letter, whose intent is just to show, through 
a simple 15 /itm-based model, what the 24 /im/15 fim 
comparison can tell us in terms of galaxy evolutionary 
properties. 

4. CONCLUSIONS 

We have discussed the importance of a comparison be- 
tween the extragalactic source counts in two different 
MIR bands for the study of the evolutionary properties of 
galaxies. We have shown what we expect in terms of vari- 
ations of the S , 24 A tm/S'i5/jm fhix density ratio with redshift 
and 24 fim flux, considering typical SEDs for the different 
populations contributing to the source counts. Then, to 
compare the observed Spitzer 24 fim source counts with 
the ISOCAM 15 fim ones, we have presented the predic- 
tion in the 24 fim Spitzer band from a phenomenological 
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evolution model based on the ISOCAM 15 /im LF of 
galaxies and AGNs. Actually, this model is the only one 
available in literature that is able to reproduce the ob- 
served 24 /im source counts without the need, at least 
in first approximation, of any a posteriori updates. We 
have also shown that the observed ISOCAM data points 
transformed from 15 /zm to 24 fim seem to agree well 
with the recently published Spitzer source counts. Our 
model suggests the appearance of a new population of 
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high redshift (z > 1.5) galaxies at 24 /im, not detected 
in the previous ISOCAM surveys, but probably the high- 
z likes of the 15 /im sources. 
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Fig. 1. — S24nm/Si5fj, m ratio as a function of redshift for the MIR populations contributing to the observed number counts: starburst 
galaxies (short-dashed line: M82 SED prototype), normal galaxies (dotted line: M51), type 2 AGN (dot-dot-dot-dashed line: Circinus) and 
type 1 AGN (long-dashed line: SED from Elvis et al. 1994). 
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Fig. 2. — Normalised differential source counts in the ISOCAM 15 fim band. As explained in the plot, data points are from several 
surveys (HDF-N, HDF-S, Marano Firback, Ultra-Deep and Deep: Elbaz et al. 1999b; Ultra-Deep lensed: Metcalfe et al. 2003; ELAIS-S1: 
Gruppioni et al. 2002; Lockman Deep and Shallow: Rodighiero et al. 04). The model curves are from Pozzi ct al. (2004) for galaxies 
(dash: starburst; dot: normal) and from Matutc ct al. (in preparation) for AGN (long-dash: type 1; dot-dot-dot-dash: type 2). 
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Fig. 3. — Normalised differential source counts at 24 fim (crosses: Deep Spitzer Surveys, Papovich et al. 2004; rilled circles: FLS, Marleau 
et al. 2004; open squares: SWIRE-N1, Shupe et al. in prep.) with the model predictions (as in figure |3 and the 15 fira data (shaded 
region) transformed to 24 (im as described in the text. 




Fig. 4. — Median values of the S24^m/ Sis,™ flux density ratio as function of 524/™, as derived by our model for all the populations 
(weighted mean: solid line; different populations: lines as in previous figures). 
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Fig. 5. — Differential rcdshift contribution to the normalized differential source counts at 24 fim. The dotted, dashed, dot-dashed, 
dot-dot-dot-dashed and long-dashed lines correspond to the 0.0-0.5, 0.5-1.0, 1.0-1.5, 1.5-2.0 and 2.0-2.5 redshift intervals respectively. 



